Mode Selection in the Spontaneous Motion of an Alcohol Droplet 
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An alcohol (pentanol) droplet exhibits spontaneous agitation on an aqueous solution, driven by a 
solutal Marangoni effect. We found that the droplet's mode of motion is controlled by its volume. A 
droplet with a volume of less than 0.1 y\ shows irregular translational motion, whereas intermediate- 
sized droplets of 0.1 — 200 /A show vectorial motion. When the volume is above 300 fA, the droplet 
splits into smaller drops. These experimental results regarding mode selection are interpreted in 
terms of the wave number selection depending on the droplet volume. 
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It is well known that an oil-water system exhibits spon- 
taneous agitation, called the Marangoni effect , which 
is driven either by a thermal gradient 0] or a chemical 
concentration gradient We focus here on the lat- 

ter system, i.e., a solutal Marangoni effect under isother- 
mal conditions. A rich variety of spontaneous agitation 
on oil- water interface has been reported both experimen- 
tallyJE 0, H S E3, EH E3, EE E3| and theoretically 
[T5L Il6j . It has been shown that the nature of the self- 
agitation is sensitively dependent on the shape and size of 
the container, i.e., on the boundary condition In 
the present study, we show that an alcohol droplet driven 
by a solutal Marangoni effect selects a certain mode of 
motion critically depending on its size. This mode selec- 
tion is interpreted in terms of competition between the 
droplet size and the critical wave number in the instabil- 
ity by the Marangoni effect. 

Different modes of spontaneous motion of alcohol 
droplets are exemplified in Figs. and where the 
droplets stained with ink (Pilot Corporation, Tokyo; 
INKSP-55-B) are agitating on 100 ml of aqueous phase 
containing 2.3 vol% pentanol in a petri dish with a di- 
ameter of 18 cm. The motion of the droplet was moni- 
tored with a high-speed video camera (RedLake MASD 
Inc., San Diego; Motion Scope PCI) at 60 frames per sec- 
ond at room temperature, and then analyzed by image- 
processing software. When the volume is less than 0.1 /A, 
the droplet maintains a circular shape during irregular 
translational motion, as shown in Fig. ^b). When the 
volume is between 0.1 and 200 fil, the droplet shows an 
asymmetric morphology and exhibits directional motion 
by maintaining its shape, as in Fig. ^c), where the speed 
is almost constant. The average speed in vectorial motion 
is greater than that in the irregular translational motion 
of the smaller droplet. When the volume is above 300 fA, 
the droplet splits into smaller droplets as in Fig. [3 Fig- 
ure |31 shows a phase diagram of the modes of droplet 
motion with a change in the volume. 

The experimental trend in volume-dependent mode- 
selection can be interpreted by considering the wave- 
length due to the Marangoni instability and the length 
scale of the droplet. When a droplet is large enough to 



accept the characteristic length of Marangoni instability, 
it exhibits morphological deformation. Once the droplet 
is deformed, the variation in curvature around its periph- 
ery causes asymmetry in the surface tension, which drives 
the directional motion of the droplet. In the experiment 
in Fig. P^ c). the concentration gradient of pentanol in 
the convex region of the droplet is higher than that in 
the concave region of the droplet due to diffusion de- 
pending on the shape of the boundary. Therefore, the 
interfacial tension gradient in the convex region is higher 
than that in the concave region and the droplet is forced 
to move toward the direction from the concave region to 
the convex region This spontaneous motion has the 
effect of maintaining the shape anisotropy. In contrast, in 
the case of a smaller droplet, the circular shape is main- 
tained under Marangoni instability, and the emergence 
of a large acceleration force is inhibited. As a result, 
the average speed of such a circular droplet is less than 
that of a larger droplet and the spontaneous agitation 
is rather irregular. The interface is significantly desta- 
bilized only when the wave number of the self-agitation 
due to the Marangoni effect is less than the inverse of a 
certain wave number, l/fc c . This mode selection is shown 
schematically in Fig. Efb). 

Next, we calculate the wave number when the inter- 
face is destabilized by a solutal Marangoni effect by con- 
sidering a small perturbation (y — Asinkx, \A\ <C 1). 
When a wave is induced along the interface, generally 
the wave tends to be damped due to interfacial tension, 
the strength of which per unit length is 71 . On the other 
hand, the concentration gradient in the convex region 
of the wave is higher than that in the concave region 
due to diffusion depending on the shape of the boundary. 
Therefore, the interfacial tension gradient in the convex 
region is higher than that in the concave region and the 
difference in the interfacial tension gradient induces wave 
growth. The force strength per unit length is 72. When 
72 > 7i, the wave becomes unstable. To simplify the 
treatment, we adopt the approximation that the shape 
of the wave consists of parts of a circle, as in Fig. EJa). 

In this case, the radius of the circle, R, and the length 
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FIG. 1: (a) Schematic illustration of the experimental setup. 
(b-1) Irregular translational motion of a droplet with a volume 
of 0.017 /il. Images of the droplet at every 1 second are shown 
on the broken line, which shows the trace of droplet motion. 
(b-2) Time trace of the speed and direction of motion, (c-1) 
(left) Vectorial motion of a droplet with a volume of 10 /il. 
Images of the droplet at every 1 second are shown on the 
broken line, which shows the trace of droplet motion, (right) 
Morphology of the droplet at 2.0 s. (c-2) Time trace of the 
speed and direction of motion. 
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FIG. 2: Spontaneous motion of a droplet with a volume of 
400 fA, exhibiting fission into small droplets. 
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FIG. 3: Phase diagram of the mode of droplet motion, 
(upper) The alcohol concentration of the aqueous phase is 
2.3 vol%. (lower) The alcohol concentration of the aqueous 
phase is 2.4 vol%. Irregular translational motion at □. Vec- 
torial motion at •. Fission into smaller droplets at o. When 
the volume was less than 0.1 yl, we hypothesize that the vol- 
ume is proportional to the 3/2 power of the cross-section of 
the droplet. In the shaded regions, the mode of the droplet 
motion is not decisive to a single specific mode. 





FIG. 4: (a) Schematic illustration of the tractable model, (b) 
Relationship between the wave number of Marangoni insta- 
bility and the size of a droplet. Closed and open circles («,o) 
correspond to the eigenvalues of the sine wave around the 
droplet, (b-1) When there are no eigenvalues where the wave 
grows, the droplet cannot break symmetry and show irregu- 
lar translational motion due to the small fluctuation along the 
periphery, (b-2) When only the fundamental wave is slightly 
unstable, the droplet exhibits vectorial motion by keeping an 
asymmetric morphology, as in Fig. 0]c-l). (b-3) When the 
wave number of the maximum instability (o) grows on a large 
enough droplet, fission into smaller droplets occurs. 



of the interface per wave length, L, are derived as 
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Supposing that 71 is proportional to d [L/ (27r/fc)] /OA, 

(3) 



where a is a constant. 

72 is proportional to the difference between the interfa- 
cial tension gradients in the convex and concave regions. 
Assuming that the interfacial tension gradients are pro- 
portional to the concentration gradients of pentanol, 
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where dc/dr\ conv is the concentration gradient perpen- 
dicular to the interface in the convex region, dc/dr\ conc 
is that in the concave region, and b is a constant. The 
dynamics of the concentration of pentanol on the surface 
of the aqueous phase, c, are given by 
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FIG. 5: Schematic representation on the interfacial ten- 
sion difference. As the interface is assumed to be cir- 
cular with the radius of R, the concentration profile is 
(co — r) Ka(\/ (a + /3) / Dr) / Ko(R) at the convex region (up- 
per), and (co — F) Io(\/ (a + (3) / Dr)/ Io(R) at the concave re- 
gion (lower). The interfacial tension is proportional to the 
gradient of the concentration at the interface; see text. 



where D is the diffusion coefficient of pentanol on the 
surface, a is the evaporation rate of pentanol, (3 is the 
dissolution rate of pentanol, c a is the concentration at 
the gas-liquid equilibrium and c w is the concentration at 
the surface-bulk equilibrium in the aqueous phase. We 
impose the Dirichlet condition at the pentanol-water in- 
terface, i.e., c = Co at the interface, where Cq is the con- 
centration of the water-rich part under the coexistence 
condition. The stationary solutions with rotational sym- 
metry of Eq. (Fig. [SJ are adopted in order to evaluate 
72. Thus, at the convex region, 
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and at the concave region, 



f dc 




\dr 


conc^ 



(co - r) 



a + P 



h \R 



D 



h ( rY 



(7) 



where V = (ac a + (3c w ) / (a + /?), R = yj(a + (3) /DR 
and Iq, I\, Kq, K\ are modified Bessel functions. There- 
fore, 
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Here, we assume < kA = i<l. When fc< 1, 

. , 2 86r^A: 2 
72 - 7i ~ -aAk H 5 — , 



(9) 
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and when k 3> 1, 
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This means 72—71 < when k 3> 1, and 72—71 > 
when k <C 1. Since 72 — 71 = at fc = k c , 



kc ~ 85 V £> 
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From the above discussion, it is clear that the pertur- 
bation increases when k < k c . This means that for larger 
r, k c is smaller, i.e., the critical volume at the mode 
transition increases. 

To verify this analysis, we performed experiments 
where T was changed by covering the dish and increasing 
the pentanol concentration in the aqueous phase. When 
the dish was covered and T increased, the motion of a 
droplet with a volume of 1.3 /d changed from vecto- 
rial motion to irregular translational motion (data not 
shown). When the concentration of the aqueous phase 
was 2.4 vol% and T increased, the critical volume at 
which the mode change is induced increases as shown in 
Fig. |21 These results correspond well to our theoretical 
expectation. 

In the present study, we found that mode selection in 
the spontaneous motion of an alcohol droplet is induced 
by a change in the droplet volume. This mode selection 
is achieved as the result of the competition in the length 
scale between the droplet size and the wavelength of the 
instability due to the Marangoni effect. The characteris- 
tic wavelength is a function of rates of evaporation and 
dissolution. The present results may shed light on how 
biological motors realize chemo-mechanical energy trans- 
duction under isothermal conditions with high efficiency 
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